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^ ' Abstract 

o : 

Excited Q-balls are studied by numerical simulations in the Minimal Super- 
symmetric Standard Model with supersymmetry broken by a gravity mediated 
' mechanism. It is found that there is a suppression factor of 0(10-2) in the 

rate at which a Q-ball can emit their excess energy compared to the rate set by 
Qh! the dynamical scale of the field, m. Furthermore, it is noted that a Q-ball can 

withstand a large amount of excess energy without losing a significant amount 
of its charge. The cosmological importance of these properties are considered 
for Q-balls in the thermal bath of the early universe. 



1 Introduction 

A scalar field theory with a global U{1) symmetry can contain a class of non- 
topological solitons, Q-balls A Q-ball is the ground state of the theory in the 
sector of fixed charge so that its stability and existence are due to the conser- 
vation of the U{1) charge. In realistic theories, Q-balls are generally present in 
supersymmetric extensions of the Standard Model |^]. In particular the scalar 
potential of the Minimal Supersymmetric Standard Model (MSSM) has several 
flat directions |p that can support lepton and/or baryon number carrying Q-ball 
solutions in its spectrum. 

In cosmology, Q-balls can play a significant role in the phenomenology of 
the early universe. The Affleck-Dine condensate is unstable with respect to 
fragmentation into Q-balls so that the small perturbations naturally present 
in the AD-condensate lead to the formation of Q-balls |Q, ^. The fate of the 
Q-balls then depends on the supersymmetry (SUSY) breaking mechanism: if 
SUSY is broken by a gauge mediated mechanism they can be absolutely stable 
and contribute to the dark matter content of the universe |^], whereas if the 
SUSY breaking mechanism is gravity mediated, the Q-balls will be short-lived 
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but can still have interesting cosmological consequences, e.g. they may explain 
the ratio of baryons to dark matter in the universe Q. The charge spectrum 
and the number of balls is dependent on the form of the potential and initial 
conditions of the AD-condensate ^, ^ . 

After the Q-balls have formed, they will be subjected to the thermal bath 
of the early universe. The effect of the background clearly depends on how 
energetic the background radiation is and how strongly the Q-ball field couples 
to it. The temperature of the universe is hence significant in determining the 
importance of the thermal background for Q-balls. 

Thermal effects on Q-balls have been previously considered via two different 
mechanisms, dissociation |5[ and dissolution |jll[|. Dissociation of Q-balls 
proceeds through the collisions of thermal particles with the "hard-core" of a 
Q-ball so that excess energy is transferred to the Q-ball. This excess energy can 
then overcome the binding energy and dissociate the Q-ball. Dissolution occurs 
due to the thermalization of the "soft edge" of the Q-ball leading to charge 
transport from the Q-ball to the background. In this work, dissolution and 
dissociation have been considered in the gravity mediated scenario by studying 
the time development of Q-balls with initial excess charge. The gauge mediated 
case has been studied by thermodynamical means in 

In addition to thermal erosion, Q-balls are also subjected to evaporation ||13[| , 
where Q-balls lose their charge via pair production at their surface. The loss of 
charge through evaporation is generally much slower than that due to thermal 
processes so evaporation becomes effective only after the universe has cooled 
enough for the thermal effects to be insignificant. 

The scalar potential in the case of MSSM with SUSY broken by a gravity 
mediated scenario is given by (omitting the ^-terms) |TT| 

= m^m\i + i^iog(^)) + (1) 

where d is the dimension of the non-renormalizable term in the superpotential, 
m is the mass scale given by the SUSY breaking scale, M is a large mass scale 
and K a negative constant of order one. The interesting cases for Q-balls are 
the d=4 u^u'^d^e'^ and d=6 u'^d^d'^ directions in the superpotential. Q-balls in 
this potential are well represented by a Gaussian approximation |jll[, where the 
Q-ball field is approximated by 4>{t, x) = e^'^^(f>oe~^ 1^ . Here lo is the so called 
Q-ball frequency and ~ can be seen from the form of the 

Gaussian approximation, Q-balls in the gravity mediated case have thick walls. 
Their charge-energy relation is well approximated by ~ mQ. 

1.1 Dissociation and Dissolution 

A particle with energy T can penetrate a Q-ball up to a distance where its 
effective mass g4>{r) is equal to its energy i.e. T ~ g4>{^)j where (/>(r) is the 
Q-ball scalar field and g the coupling of the thermal particle to the ball. The 
thermal particle entering the ball will then lose some of its kinetic energy to 
the classical Q-ball field and come to a temporary halt. The excess energy 
transferred to the Q-ball can then be radiated as scalar field waves or lumps of 
charge, or some combination of these. 

The rate of the number of particles hitting a Q-ball of radius i? in a thermal 



path of temperature T is 



^ = ~^A.R'T\ (2) 
dt vr"^ 

where g{T) is the effective number of thermal degrees of freedom couphng to 
the Q-bah scalar. Typically we can set g{T) ~ 10^. 

The question of how much of the energy of an incoming particle will be 



transferred to the Q-ball was considered in [11|. Assuming that the Q-ball field 
can reconfigure itself over a time scale of m~^, the Q-ball can typically absorb 
most of the energy of the incoming particle. Hence, in each collision an energy 
of ^tT, where 7t ^ li is transferred to the Q-ball. The rate at which energy is 
being added to the Q-ball by the thermal background is therefore given by 

^=,,T^ = m4,,H=r'. (3) 

dt dt TT 

To avoid dissociation, a Q-ball should be able to radiate the excess energy 
as scalar field waves at a rate faster than given by (^). Otherwise, energy will 
be accumulated in the Q-ball until it begins to dissociate. 

In addition to dissociation, a Q-ball can lose its charge by dissolution where 
charge is removed from the "soft edge" of the Q-ball. For T > m, the thermal 
background quarks have a much smaller mean free path, Xmfp ~ kg/T {kg ^ 6) 
than the width of the soft edge of a Gaussian Q-ball, 5r R/{2(5), where j3 = 
ln(g'(/)o/3T)^/^. The charge contained in the soft edge is approximately Qsoft ~ 
A-KUj{T/gf{l3Rf6r Q. The background quarks can therefore thermalize the 
soft edge of a Q-ball at high temperatures. The time taken for a charge to leave 



the soft edge by a random walk process is |11] 



^'^'^^W k'q Al3'^m?\K\kq' ^ ' 

The rate at which charge can be removed then depends on how quickly the Q- 
ball can reconfigure itself and replenish charge in the soft edge. If the charge is 
replenished quickly enough, the rate of charge loss is dQ/dt ^ Qsoft/Tdi where 
Qsoft is the charge contained in the soft edge. However, if the Q-ball does not 
reconfigure and hence replenish the charge quickly enough, the rate of charge 
loss will be suppressed. 



2 Numerical simulations 
2.1 Preliminaries 

To study the dynamics of excited Q-balls, they have been simulated numerically 
on a 2+1 dimensional lattice. Obviously, a 3-1-1 dimensional simulation would 
be more realistic but it is expected that a 2-1-1 dimensional simulation gives a 
fair approximation to the more realistic case while greatly reducing the required 
CPU time. 

The equation of motion in flat space corresponding to the potential (||) is 
* + V^^ + + ^ + A',og(ffl!)l + «^|$|^(-^)* = 0. (5) 



(We can neglect the expansion of the universe for our purposes since we are 
interested in the behaviour of a single, isolated Q-ball.) 

Since our goal is to address questions relevant to Q-balls subjected to a ther- 
mal background radiation, it is necessary to consider exactly what to simulate 
while taking into account the numerical and practical problems that may arise. 

Simulating the collision of a single particle of energy T with the Q-ball is 
cumbersome due to the vastly different energy scales involved. In a realistic 
cosmological setting, the energy of a Q-ball is typically much larger than T. 
Keeping track of such a tiny change in energy would require a simulation where 
the Q-ball field is let to evolve at very high accuracy while keeping the time 
and spatial steps very small. This in order requires a large lattice and several 
orders of magnitude more time steps. Also simulating a single particle hitting 
the Q-ball is not realistic since a Q-ball in the thermal bath of the early universe 
is subjected to several particles interacting with it simultaneously. 

Another approach would be to fully simulate both a generic background field 
and the Q-ball field. The difficulty in this approach is to keep the background 
field at a constant temperature at large distances i.e. to keep the system from 
reaching a global equilibrium state due to the smallness of the lattice compared 
to the Q-ball size. By using a much larger lattice this approach would probably 
be more feasible but would again require more CPU time. Hence an approxi- 
mate method has been adopted: instead of simulating the whole thermal erosion 
process, we simply study the evolution of a Q-ball with excess energy added in 
the form of a random perturbation, (/> = (/>o + 6(j), where 0o is the Q-ball field and 
6(p is a stochastic term. Such a perturbation adds to the energy of the Q-ball 
while keeping the contribution to the charge small. The exact form of the per- 
turbation is not important; the relaxation times are insensitive of how the extra 
energy is added to the Q-ball. For example, if energy is added asymmetrically, 
it has been noted that the fiuctuations of the field quickly distribute around the 
Q-ball before excess energy begins to radiate. 

After the perturbation is added to the initial configuration, the field is al- 
lowed to develop until the Q-ball becomes stable again i.e. excess energy and 
possibly charge is radiated from the ball. To avoid any excess energy and pos- 
sible charge radiated by the Q-ball from entering the system due to continuous 
boundary conditions, the lattice has dissipative zones near the boundaries (to 
minimize back-scattering from the boundaries of the dissipative regions the dis- 
sipation term is increased gradually). For the numerical simulations the field 
was decomposed into real and imaginary parts, $ = -^{<Pi + ^^2), and re-scaled 

along with the spatial and temporal unit, (fi = r = mt, ^ = mx. The 

parameter values wc used in the numerical simulations were m = 100 GeV, K = 
-0.1, M = 10^ GeV. The calculations were done mainly for the d = 4 case but 
we have checked that the d = 6 case does not significantly differ from these 
results. We have also verified that varying K in the typical range —(0.01 — 0.1) 
docs not markedly affect the relaxation rate. The lattice size was typically 
200 X 200 and the lattice parameters were chosen as = 0.6, At = 0.3 (a 
number of simulations with varying values of and Ar were done to ensure 
that the exact values of the lattice parameters do not affect the results). 



2.2 Numerical results 



We have studied numerically the relaxation of Q-balls for different Q-ball charges 
and excess energies. For each set of parameter values a number of simulations 
were run to account for the effect induced by the stochastic term i.e. a different 
seed for the pseudo random number generator in the stochastic term was used. 

A typical relaxation process is shown in Fig. || where the amplitude of the 
scalar field is plotted. The initial configuration, with to = 0.80m corresponding 
to Q = 9100 and E = 7600 (in lattice units), is shown in Fig. |^(a), is put on 
a 200 X 200 lattice and the field is let to evolve. Note that the images are not 
taken at equal time intervals: the time units at which snapshots are taken are 
Ta = 0, Tf, = 20, Tc = 40, Td = 990. From the figures it can be seen how the 
Q-ball quickly radiates the excess energy in the form of scalar field waves which 
propagate spherically from the ball. The waves are dissipated at the dissipation 
zones as they reach the boundaries of the lattice. Radiation gradually decreases 
and in the last panel one can see virtually no radiation. Prom Fig. |l|(d) one 
can also note how the Q-ball has slightly moved from its initial position in the 
center of the lattice due to a spatial asymmetry present in the initial conditions. 

The energy and charge of the Q-ball in lattice units have been plotted in Fig. 
|2|. From the figure it can be seen how energy is emitted from the ball while the 
charge remains virtually constant, even though a significant amount of excess 
charge is introduced to the Q-ball initially. The slight fluctuation of charge at 
the beginning of the simulation is shown in the close up in Fig. §(b). 

To study how the energy is emitted from the Q-ball a number of different 
functions has been fitted to the energy curves. It is found that the energy curve 
is well described by 

E{t) = A[e{To - r) + e{T - ro)e-^("-"«)] + ^o, (6) 

where A, A, Eq are the fitting parameters and 9(t) is the step function. The 
presence of a non-zero TQ-term implies that the excess energy first spreads around 
the ball before it is emitted. 

The rate at which excess energy is emitted has been studied for iv = 0.78m — 
0.98m, which in two spatial dimensions corresponds to Q-ball charges in the 
range Q ~ 10^^ — 10^*^ (note that smaller values of uj correspond to larger 
charges). The values of A have been plotted in Fig. |3| as a function of the 
unperturbed Q-ball energy (in lattice units) . For each value of oj the maximum 
amplitude of the perturbation 6(j) was varied, \6(j)\ < a(f){0), where a was varied 
in the range 0.01 — 0.10. The maximum value of a for each to was chosen so that 
the charge of the Q-ball varied less than 1% from its initial value during the 
relaxation process. Ten simulations were run with a different initial stochastic 
perturbation for each a. 

The amount of excess charge is obviously a function of a and hence it varies 
greatly for each value of lo. For example, for the lo = 0.78m case, the values of A 
have been calculated for increases in energy between 1% and 40% of the energy 
of the ground state. However, as from Fig. |^ can be seen, the variations in A are 
small compared to the changes in energy. The values of A increase slightly with 
smaller energy so that smaller Q-balls radiate their excess energy more quickly 
than larger Q-balls. 

We have also studied the effect of adding larger amounts of excess energy to 
a Q-ball. The time development of energy and charge of a Q-ball with perturba- 
tions of different magnitude for to = 0.80m is shown in Fig. ^ The ground state 



has charge Qq = 9100 and energy Eq = 7600. From the figure it can be seen how 
energy still decreases exponentially even with very large initial excess energy. As 
the amount of excess energy grows, the Q-ball emits more charge along with the 
energy but still retains most of its charge in a single Q-ball. In the case shown 
in Fig. ^, more than 75% of the initial charge remains in the Q-ball even if the 
excess energy is ~ 10^. A Q-ball can hence survive fragmentation even with 
initial excess energies of the order of ^ 10 times its initial energy. 



3 Discussion 

From the numerical results we have seen that the energy loss experienced by an 
excited Q-ball is well described by an exponential function, EQ{t) = —\{EQ{t) — 
Eqq ) , where Eq^ is the ground state of the Q-ball. 

To study the dissociation of Q-balls in the early universe we need to consider 
the energy increase due to the thermal background. Since the rate at which 
energy is being added to the Q-ball due to the thermal background is given by 
(|3|), we can model the energy of the Q-ball in the early universe by 

EQit) = -XiEgit) - Eq,) + ^47Ti?'r^ (7) 

Here we can approximate that R is approximately independent of the charge (and 
energy) of the Q-ball, assuming that the energy increase is not large enough to 
dissociate the Q-ball. 

The temperature, T, decreases with time due to the expansion the universe 
so that in a radiation dominated universe T^t ~ Mpi and in a matter dominated 
universe T^/^t ~ Afpi- In general one can therefore write (|^ as 

EQ{t) = -XiEQ{t)-EQ,) + ^, (8) 

where j4 is a constant and B = 2 (radiation dominated) or S = | (matter 
dominated). Integrating once with respect to time one gets 

dx (9) 

to 

where Iq is now the time at which the Q-ball is formed. 

The relevant features of the time evolution of Q-ball energy can be extracted 
from (^). As t ^ CX3, the energy of a Q-ball tends to the ground state (assuming 
that the Q-ball does not dissociate during the process) i.e. EQ{t oo) = 
Eqq. If there exists an extremum tc, E{tc) = 0, any subsequent extremum 
points, E{t) = 0, will occur at a lower energy i.e. E{tc) > E{t) for B > 1. 
Hence the value of E{tQ) determines the features of the time evolution of a Q- 
ball: If i?Q(to) < ^/{XtQ) + Eqq, E{to) is positive and the Q-ball energy will 
initially increase until the energy loss from the Q-ball exceeds the contribution 
from the background. The increase in energy is clearly bounded from above by 
Atl~^/{B - 1). If, on the other hand, S(to) < 0, the energy of a Q-ball will 
monotonously decrease with time. 

If one assumes a situation where Q-balls are formed before reheating due to 
the inflaton decays and are approximately at their ground state, EQ{t = to) ~ 



Eqq, the energy of a Q-ball will initially increase, E{t) > 0. The temperature of 



the universe after reheating can be approximated by T{t) = y t > tRH = 
Mpi/Tp^jj, where Trh is the reheating temperature and tf/jj is the reheating 
time. The maximum energy of a Q-ball clearly depends on the exact value of 
parameters. However, for typical parameter values, numerical work shows that 
for Trh ~ 10^ - 10^ GeV the increase in energy is well approximated by 

> tRH) « EQitnn) + ^A^tR'T^h, (10) 

which corresponds to the point at which the rates at which energy is emitted 
and absorbed are equal. Clearly, ([lO| ) is valid only if the Q-balls emit their 
excess charge quickly enough compared to the rate of energy increase due to the 
energetic background. 

From the simulations we know that to significantly remove charge from a 
Q-ball we need to add a perturbation whose energy is comparable to the energy 
of the Q-ball itself. Therefore, to avoid significant loss of charge the energy of a 
Q-ball is bounded from below by 

EQ^mQ> ^4jtR^T^h, (H) 
vrAa 

where a is the fraction of excess energy that is required to significantly remove 
charge from the Q-ball. Using R ~ m~^\K\~^/'^ , K = —0.1, 7^ ~ 1 {i-e. all 
of the energy of an incoming particle is transferred to the Q-ball), g{T) ~ 100, 
A ~ 0.02m and a ~ 1 we get 

Tfl// < 6 X GeV. (12) 

On the other hand, if Q-balls are still at an excited state at the time of 
reheating and the reheating temperature is low, the rate of energy loss from a 
Q-ball can be greater than the rate of energy increase due to reheating. The 
energy of a Q-ball can then be monotonously decreasing with time during the 
reheating process and dissociation due to the energetic background will not play 
a major role in the evolution of Q-balls. 

In addition to dissociation, dissolution can also be a significant factor in 
determining the fate of Q-balls in the thermal bath of the early universe. For 
dissolution to be important the rate of charge loss from the soft edge, r^"^, 
must be smaller than the rate at which charge is replenished. The rate which a 
Q-ball can reconfigure itself and replenish charge in the soft edge now becomes 
important. The numerical simulations show that the Q-ball emits excess charge 
at a slow rate, F ~ Am compared to dynamical scale of the scalar field, m. One 
can then speculate that there is a suppression factor in how quickly charge is 
replenished in the soft edge, ~ (^m)^^. 

For efficient dissolution we need t^^ ^ ''"r"^ so that using (Q) we get a lower 
bound for the background temperature for , 

T>4.l3'^m\K\kqK-^. (13) 

Using typical values, we get T > I X 10^ GeV so that the presence of a suppres- 
sion factor can significantly increase this bound. For example, assuming that 
K ~ A, we get T > 10^ GeV. 



At temperatures high enough for efficient dissolution, charge in the soft edge 
is replenished quickly enough and dissolution is efficient. The rate of charge loss 
is then given by, 

^ ^ ^ = 8^co(3^k,Rg-'T, (14) 
at Td 

where @ has been used. Using T = {M-pi/lkxtY^"^ and integrating, the total 
loss of charge from the time of reheating until temperature has decreased to T, 
is given by 

1 1 

The temperature at which the loss of charge due to dissolution is no longer 



AQ = 8TTu;f3^kgRg-^kj,^Mpi{- - - — ). (15) 



efficient is given by (|l3| ) and hence (15) is bounded from above by 



~ g'^kT'm'^\K\'^/'^ 

For typical values to ^ m, /3 ~ 4, g l,kT 17 and K = —0.1 we get 
AQ ~ 0.5K(Mpi/GeV). One should also note that the rate at which charge 
leaves the soft edge of a Q-ball by random walk will in reality be less than Td~^ 
since random walk transports charges in all spatial directions and the rate at 
which charge drifts away from the Q-ball will depend on the gradient of the 
charge density. Furthermore, charge diffusion away from the edge of the Q-ball 
needs to be efficient to prevent charges being transported back into the Q-ball. 
These considerations can further suppress the dissolution process. 



4 Conclusions 

We have studied the effect of adding excess energy to a Q-ball by numerical 
means. Prom doing a set of simulations we have found that the excess energy 
is lost as E{t) ~ e~^'^ ~ e"^*"*, where A ~ 0.02 — 0.06 in the studied parameter 
range. From naive arguments one would expect that the relaxation rate of 
the Q-ball would be of the order of m. The numerical simulations show that 
there is a suppression factor 0{W~^) in how quickly Q-balls radiate their excess 
energy. We expect that this suppression factor would also be present in full 
3D simulations even though the Q-balls are likely to radiate their excess energy 
slightly more quickly due to the presence of an additional degree of freedom. 
We have also observed the slow relaxation rate of excited Q-balls when studying 
Q-ball collisions [||] . 

In addition to the fact that Q-balls tend to radiate their excess energy more 
slowly than what one might expect, it has also been found that the Q-balls are 
surprisingly robust. Even when the excess energy added is a significant fraction 
of the energy of the unperturbed Q-ball and clearly exceeds the binding energy, 
niQ — E, only a very small fraction of the charge is lost in emitting the excess 
energy. Only when the excess energy is much larger than the Q-ball energy we 
start to see a significant loss of charge. One can gain an insight into this by 
looking at the effective potential of a single "free" (j) particle inside the Q-ball: 
the particle is typically in a deep potential well with a high barrier that also 
suppresses quantum tunneling. 

On the basis of the small rate of radiation of the excess energy, one can 
speculate that also the reconfiguration rate may be suppressed. This would 



then reduce the rate at which a Q-ball can absorb energy from the thermal 
background and suppress the charge loss by dissolution. Especially the effect 
on dissolution is significant since a suppression factor in the reconfiguration rate 



affects linearly the lower bound for efficient dissolution (13) and hence also the 
upper bound for charge loss ([l6|) . 

The slow relaxation rate and the robustness of Q-balls have an effect on how 
well they can avoid thermal erosion in the early universe. The slow relaxation 
rate allows for the excess energy to build up in a Q-ball but on the other hand the 
excess energy required to significantly reduce the charge of Q-ball by dissociation 
is large. In a realistic cosmological setting one should also take into account the 
effect of preheating due to infiaton decays and the decays of the remains of 



the AD-condensate. However, on the basis of the arguments presented in 
preheating effects should not significantly affect the bounds for the Q-ball charge 
and reheat temperature. 

Q-balls formed from the Affleck-Dine condensate can survive the thermal 
background of the early universe if the reheat temperature is not too large. 
On the basis of the results presented in this work Q-balls with charges ~ 10^'^ 
should survive reheat temperatures around Trh ~ 10^ GeV. The question of 
how quickly a Q-ball can absorb energy and replenish charge in its soft edge can, 
however, change theses limits so further work is still needed to understand the 
exact dynamics of how Q-balls erode thermally. 
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Figure 1: A typical Q-ball relaxation process 
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Figure 2: Q-ball energy (a) and charge (b) as a function of r, a; = 0.80m 
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Figure 3: A for a; = 0.78m — 0.98m as a function of Q-ball energy 
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Figure 4: Q-ball energy (a) and charge (b), a; = 0.80m, Qq = 9100, Eq = 7600 



